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In this paper we review the latest achievements of density functional theory in understanding 
the physics of diluted magnetic semiconductors. We focus on transition metal doped III-V semi- 
conductors, which show spontaneous ferromagnetic order at relatively high temperature and good 
structural compatibility with existing III-V devices. We show that density functional theory is a 
very powerful tool for i) studying the effects of local doping defects and disorder on the magnetic 
properties of these materials, ii) predicting properties of new materials and iii) providing parame- 
ters, often not accessible from experiments, for use in model Hamiltonian calculations. Such studies 
are facilitated by recent advances in numerical implementations of density functional theory, which 
make the study of systems with a very large number of atoms possible. 
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I. INTRODUCTION 



The possibility of using the spin degree of freedom as 
well as the electronic charge for electronic applications 
( "spintronics" ) has received great attention in the last 
few years [Q. A prototypical application of this con- 
cept is the giant magnetoresistance effect (GMR) 
where the electrical resistance of a magnetic multilayer 
is changed by the application of a magnetic field. GMR 
sensors are already commercially available for read heads 
in high density data storage devices, with performances 
increased by one order of magnitude with respect to con- 
ventional sensors. However current GMR sensors do not 
exploit completely the potential of spintronics. In exist- 
ing electronic devices (for example personal computers) 
there are two main elements: the logic components and 
the data storage device. The former are transistors based 
on semiconductor technology, while the latter is essen- 
tially a metallic magnetic film. Of course the ability to 
combine both logic elements and data storage in the same 
device will open completely new possibilities, with huge 
potential applications. Hybrid structures, where mag- 
netic metals are used to inject spin electrons into semi- 
conductors have been shown to be problematic ^, 1). 
This is due to the large mismatch between the resistivities 
of semiconductors and metals, which seriously precludes 
effective spin- injection [||. Although this fundamental 
obstacle can be overcome, for instance by large contact 
resistances [fTf, it is natural to turn attention towards 
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other directions and try to create all-semiconductor de- 
vices. These, of course, need magnetic semiconductors. 

At present several types of diluted magnetic semicon- 
ductors (DMSs) are available |^. Among these, the 
DMSs based on III-V semiconductors are particularly 
important because of their compatibility with existing 
Ill-V-based technology. III-V DMS are obtained by low- 
temperature molecular beam epitaxy (MBE) deposition 
of III-V semiconductors with a transition metal such as 
Mn. The non-equilibrium growth is necessary to pre- 
vent the formation of additional phases and in general 
only low concentrations of transition metal ions can be 
incorporated in the non-magnetic matrix. Nevertheless, 
despite the low concentrations, the systems develop long 
range ferromagnetic order with remarkably high Curie 
temperatures, Tc- For the known Ill-V-based DMS the 
highest Curie temperatures obtained are: Tc ^ 30K for 
(In,Mn)As |, §, Tc^UOK for (Ga,Mn)As Hand, very 
recently, a report of T^^MOK for (Ga,Mn)N |[lO|. More- 
over, although this paper is focused only on the III-V 
DMS, it is worth mentioning that ferromagnetism with Tc 
often above room temperature has been found in several 
other DMSs, including Gei_a;Mn^ Cdi_^Mn3;GeP2 
Tii_,Co,02 |l| and Zni_,Co,0 @. 

(Ga,Mn)As, although it has a Tc far below room tem- 
perature, is the most widely studied DMS. This is due to 
its structural compatibility with most epitaxially grown 
III-Vs, which makes it ideal for building heterostructures 
and exploring new device concepts. Examples that have 
been already achieved include spin injection into het- 
erostructures and field-effect control of the ferromag- 
netism 0. 

Three important features underlie the ferromagnetic 



2 



order of (Ga,Mn)As, and these are shared by the other 
Mn-doped III-V DMSs: i) Mn'^+ ions substitute the for 
Ga^+ cations in the zincblende lattice providing local- 
ized spins (5* = 5/2 in Gai_a;Mn2;As and Ini_a;Mna;As), 
ii) there are free holes in the system although the actual 
concentration is much smaller than the Mn concentra- 
tion (despite the nominal valence suggests that the two 
concentrations should be identical), iii) the hole spins 
couple antiferromagnetically with the Mn spins, due to a 
dynamic p-d coupling. Then the ferromagnetic behavior 
can be reasonably well described by the Zener model , 
in which antiferromagnetic exchange coupling partially 
spin polarizes the holes, which in turn cause an align- 
ment of the local Mn spins. Within the Zener model, the 
interaction Hamiltonian between the hole spin s and the 
Mn spin S is 
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where Nq is the concentration of the cation sites and /3 
is the p-d exchange integral. Nq(3 is usually called the 
exchange constant. If one simply uses the mean field 
approximation |T^, in which the magnetizations of 
both carriers and Mn ions are considered to be uniform 
in space, we obtain an expression for Tc 
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where Xs is the magnetic susceptibility of the free car- 
riers (holes in this case), g* is their g- factor, ks is the 
Boltzmann constant and /is the Bohr magneton. 

The result of equation can be greatly refined by 
including a detailed description of the band structure of 
the underlying non-magnetic semiconductors or by incor- 
porating correlation effects going beyond the mean field 
approximation. However, it is important to stress that 
the use of these models always requires parameters that 
are often difficult to extract from the experiments. For 
example the experimental value of the exchange constant 
NoP vary in the range 1-3 eV Q |o[ |2|, ||, ||. 

In addition, models based on any kind of mean field 
approximation naturally fail in describing local effects, 
which occur on the atomic scale. The ferromagnetism of 
(Ga,Mn)As is very sensitive to the sample "history" , such 
as the growth conditions and eventual after-growth 
processing [^s], ^ . Since the growth dynamics certainly 
affects the microscopic configuration of the samples, this 
suggests that knowledge of the local chemical environ- 
ment is crucial for understanding and modeling the prop- 
erties correctly. 

These considerations show that it is essential to have a 
microscopic theory providing information to the simpler 
mean field-like models. Density functional theory (DFT) 
p7| , p8[ is to date the most efficient and accurate micro- 
scopic theory for describing the electronic, magnetic and 
structural properties of the ground state of electronic sys- 
tems with a large number of degrees of freedom. Recent 



advances in the numerical implementations, mainly con- 
centrated in the use of improved pseudopotentials [ p9| , 
of efficient basis sets and of order TV approximations 
PH , make possible the study of systems containing sev- 
eral hundreds of atoms. Such computational capabilities 
are required to study the DMSs in the low concentration 
limit. The main aim of this paper is to provide a review 
of the achievements of DFT in describing the properties 
of the DMSs. In particular we will show that DFT is an 
invaluable tool for i) studying the effects of local dop- 
ing defects and disorder on the magnetic properties, ii) 
predicting properties of new materials and iii) providing 
parameters, often not accessible from experiments, for 
use in model Hamiltonian calculations. 

The remainder of the paper is organized as follows. 
First we will briefly overview the most recent devel- 
opments in density functional theory. In the follow- 
ing section we will discuss the structural properties of 
(Ga,Mn)As and explain why zincblende MnAs cannot 
be grown. Moreover we will also look at spin-orbit ef- 
fects and explain why this introduces only minor quan- 
titative changes in the exchange constant of (Ga,Mn)As. 
Then we will move to the low dilution limit, calculating 
the exchange constants and discussing the limitations of 
mean-field models. In the remaining sections we will con- 
sider the effects of the local microscopic configuration of 
the Mn ions and possible intrinsic defects on the ferro- 
magnetism of (Ga,Mn)As. In particular we will look at 
the role of intrinsic defects and at the transport proper- 
ties of digital ferromagnetic heterostructures (DFH) . 
Finally we will overview theoretical predictions for new 
materials and then we will conclude. 



II. DENSITY FUNCTIONAL THEORY 

Since its introduction in the 1960s ^ density func- 
tional theory has evolved into a powerful tool that is 
widely used in condensed matter theory and computa- 
tional materials for the calculation of electronic, mag- 
netic and structural properties of solids. The method 
has been remarkably successful in predicting, reproduc- 
ing and/or explaining a wide variety of materials phe- 
nomena. Specific examples range from early predictions 
of p hase transitions in silicon as a function of pressure 
p3| , to determination of stable and meta-stable adsorp- 
tion geometries on metal surfaces [|4| as well as many 
successes in understanding the behavior of magnetic ma- 
terials, including those described in this work. 

The density functional formalism is based on the the- 
orem that for an interacting inhomogeneous electron gas 
in a static external potential, w(r), there exists a uni- 
versal functional of the density, F[p{r)], independent of 
v{r), such that the expression 
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has as its minimum value the correct ground state energy 
associated with v{r) pTf . 

The true density, p(r), can in principle be exactly 
obtained from the solution of an associated single- 
particle problem, whose effective single-particle poten- 
tial, Veff[p{v)], is a unique functional of p(r) As a 
consequence, the many-electron ground state reduces to 
that of a one-electron Schrodinger equation: 
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The so-called Kohn-Sham wavefunctions, <f>i, are single- 
particle eigenfunctions that are strictly meaningful only 
for determining p(r), and the Kohn-Sham eigenvalues. 
Si, are the derivatives of the total energy with respect to 
the occupation of state i. Note that they are not strictly 
related to single particle excitation energies, although the 
Kohn-Sham band structure can sometimes be a useful 
tool in the interpretation of photoemission (or similar) 
data. 

For an arbitrary density there is no simple exact ex- 
pression for the exchange-correlation energy, E^c, and 
so to make progress, the so-called local density approx- 
imation (LDA) is often made. Within the LDA, E^c is 
written as: 



Exc[p\ = j p(r)£ 
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where Cxc is the exchange-correlation energy per electron 
of a uniform interacting electron gas of the same density, 
p. The LDA is strictly valid only if p(r) is slowly varying, 
and many extensions exist which give improved accuracy 
for systems with localized electrons. 

Density functional calculations for magnetic materi- 
als became widespread in the late 1970s, with a num- 
ber of studies of third and fourth row transition metals 
p5| , |3^ , |37| | . These studies established that the local den- 
sity approximation gives results that are in reasonable 
agreement with experiment for quantities such as cohe- 
sive energy, bulk modulus and magnetic moments, pro- 
vided that spin polarization is included explicitly, by ex- 
tending the LDA to the local spin density approximation 
(LSDA). They also noted, however, that the calculated 
properties are very sensitive to details of the structure 
and magnetic ordering, which can lead to discrepancies 
between the LSDA results and experiment. The most 
notorious of these is the well known prediction of the 
incorrect ground state of iron (face centered cubic and 
antiferromagnetic, rather than the correct body centered 
cubic and ferromagnetic) by the LSDA. 

A number of technical developments have facilitated 
the study of magnetic materials, perhaps the most im- 
portant being the introduction of the fixed spin moment 



(FSM) method |||, |3|l. In the FSM method the ground 
state of a constrained system with a fixed magnetic mo- 
ment is calculated. Not only this does speed convergence, 
but the total energy surface in magnetic moment/volume 
space can be determined, giving additional information 
particularly about metastable magnetic phases. Also, 
implementation of Gaussian smearing [Q and related 
schemes have helped to speed convergence of calculations 
for magnetic metals with partially filled d bands and com- 
plex Fermi surfaces, in which it is difficult to carry out 
integrals over the occupied part of the Brillouin Zone. 

In parallel with these technical developments, exten- 
sions and improvements to the LSDA have also been 
explored. The usual generalized gradient (GGA) and 
weighted density (WDA) approximations that give im- 
proved results for non-magnetic systems do not tend to 
give systematic improvement for magnetic materials, al- 
though the GGA does at least predict the correct ground 
state for iron. For more information about these approx- 
imations see reference and references therein. Meth- 
ods such as the LDA+U and self- interaction- 
correction iQ are specifically tailored to treat strongly 
correlated systems, and therefore are more appropriate 
for magnetic systems with narrow d or f bands. 



A. Codes available 

There are many excellent computer programs available 
today for performing density functional theory calcula- 
tions. These use a range of different methodologies, have 
different specialties, and are widely varying in cost (both 
in dollar amount and in their computer requirements). 
In this Section we describe some of the most popular 
programs. 

Many DFT programs are based on the so-called plane 
wave pseudopotential (PWPP) method |3^, in which the 
wavefunctions are expanded in a plane- wave basis, and 
the electron-ion interaction is modeled by a pseudopo- 
tential. Plane wave basis sets offer many advantages 
in DFT calculations for solids, including completeness, 
an unbiased representation, and arbitrarily good conver- 
gence accuracy. Publically available PWPP codes in- 
clude DoD Planewave [^s) which is a general purpose 
scalable planewave basis density functional code that 
treats insulators, semiconductors, metals and magnetic 
materials, ABINIT [|6), which allows both DFT and 
density functional perturbation theory calculations, and 
Spinor which extends the usual LSDA formalism to 
include spin-orbit coupling and generalized non-coUinear 
magnetism. Both ABINIT and Spinor are PWPP 
codes and are available under the GNU General Public 
License |Q . One of the most popular "semi-commercial" 
PWPP codes is the VASP package [|9) developed at the 
University of Vienna, Austria. VASP allows DFT and 
molecular dynamics calculations, and is quite fast be- 
cause of its use of ultra-soft pseudopotentials. The de- 
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velopers charge a nominal fee for the source code, and re- 
quire authorship on the first pubhcation using the code. 
There are also a number of fully commercialized density 
functional codes that are targeted in large part at chemi- 
cal and pharmaceutical companies. For example accelrys 
markets the PWPP CASTEP code ||]. 

For systems involving a large number of atoms in the 
unit cell, plane-wave based DFT algorithms are not ideal 
because of the large computational overheads involved. 
For this purpose it is convenient to use codes based on 
localized atomic orbital basis sets, although their numer- 
ical implementation is usually quite complicated. Most 
of the results of this paper are obtained with the code 
SIESTA ^2, which combines pseudopotential tech- 
niques with a pseudoatomic orbital basis set |3^. The 
code is highly optimized to deal with large systems, and 
both efficient order N methods and molecular dynam- 
ics tools are available. The developers require an initial 
collaboration and co-authorship on works produced by 

SIESTA. 

Traditionally however, magnetic materials have been 
studied using all-electron methods with mixed basis sets, 
such as the linear augmented plane wave (LAPW) 
linear muffin tin orbital (LMTO) jsj] or Korringa-Kohn- 
Rostoker (KKR) ||5^ approaches. Again many codes 
are available. For example, the Vienna University of 
Technology produces a full-potential LAPW code called 
WIEN97 which is a highly accurate, all-electron 
code that includes relativistic effects. A small ($350 at 
press time) fee is charged. The Stuttgart LMTO program 
^Tj is a fast and efficient tool for calculation of charge- 
and spin-self-consistent band structures, partial densities 
of states, Fermi surfaces, total energies, and the partial 
pressures. In addition, the program delivers tools for 
analyzing the electronic structure and chemical bonding 
such as orbital-projected band structures, crystal orbital 
Hamiltonian populations and electron localization func- 
tions. 



III. ZINCBLENDE MnAs 

As we mentioned in the introduction, so far Mn has 
been incorporated in GaAs only at concentrations smaller 
than 7%. Considering that the Zener model predicts that 
Tc increases with x (see equation (^), this poses a severe 
limit to the highest Curie temperature obtainable. It is 
therefore of great interest to investigate if there are some 
conditions, which allow the growth of (Ga,Mn)As with 
higher Mn concentration or even ultimately zincblende 
MnAs. In this Section we describe DFT calculations 
from the literature that investigate both the stability and 
the electronic and magnetic properties of the zincblende 
MnAs. 

Several groups have calculated the electronic structure 
of zincblende MnAs using both the LSDA g, ||, |0|, |6|] 
and the GGA 162] approximations for the exchange cor- 



relation potential. The calculated lattice constant, oo, 
for zincblende MnAs is found to be around 5.6-5.7 A. 
This result, also confirmed by relaxation calculation of 
(Ga,Mn)As is quite appealing since a hypothetical 
zincblende MnAs appears to have a lattice constant very 
close to that of GaAs (aq = 5.6533A). These results seem 
to conflict with the extrapolation to x=l (assuming Ve- 
gard's law) of the lattice constant of Gai-a^Mn^^As and 
Ini_2;Mn2;As 1^, which give respectively ao=5.98A and 
ao=6.0lA. However it is worth noting that the lattice 
constant has been measured only for a; ^ 1 and that a 
linear extrapolation to a; = 1 may be not valid. More- 
over in real (Ga,Mn)As samples the presence of intrinsic 
defects (mainly As antisites Asca) can play an impor- 
tant role in determining the structural properties. This 
is clearly demonstrated in low-temperature annealing ex- 
periments ll^, ^ , where the lattice constants of samples 
with the same Mn concentration, change depending on 
annealing conditions (temperature and time). 

It is also important to point out that generally the 
LDA approximation underestimates the equilibrium lat- 
tice constant, in particular if strong p-d hybridization is 
present. This is probably the situation for zincblende 
MnAs, though the good agreement with experiments for 
a different MnAs lattice structure (the NiAs-type struc- 
ture), gives confidence in the LSDA results. Furthermore 
a recent LDA-I-U calculation |Q| (LDA-I-U usually cor- 
rects the tendency to ovcrbinding of the LDA) finds a 
lattice constant for MnAs very similar to that of GaAs 
|]65| . This suggests that LDA provides a good description 
of, at least, the structural properties. Finally we must 
point out that a much larger lattice constant (ao=5.9A) 
has been found by Shirai et al. |6^ within the LSDA. 
However this is probably due to an artifact of the mini- 
mization procedure used. 

Turning our attention to the electronic properties, in 
figure we present our calculated band structure for 
MnAs at the LDA equilibrium lattice constant ao=5.7A 
from referenc e |6C[| . The results are obtained with the 
code Spinor |}47[. Calculations by other groups at sim- 
ilar lattice constants give similar band structures [|62[ . 
We first note that, excluding the presence of the Mn d 
bands, the bandstructure closely resembles that of the 
non-magnetic III-V semiconductors. If we consider the 
majority band first, we can easily identify the As p va- 
lence band (first Fis point above E-p) and the first of the 
conduction bands (first Fi point above E-p). However the 
strong interaction with the Mn d states pushes the former 
toward higher energies and they become half filled. The 
Mn d bands, which are split into the doubly degenerate 
e band (ri2) and the triply degenerate t2 band (Fis) are 
below the Fermi energy and entirely occupied. In con- 
trast, in the minority band there is a large split between 
the t2 and the e states, which gives rise to a large gap in 
the band structure. 

It is very important to note that at this lattice constant 
the Fermi energy in the minority spin band cuts through 
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FIG. 1: Band structure for zincblende MnAs at the LDA en- 
ergy minimum (ao=5.7 A). The figure on the left corresponds 
to the majority spin and the one on the right to the minority. 
The horizontal line denotes the position of the Fermi energy, 
which has been chosen to be eV. 



the edge of the Mn-c? (e) states. This suggests that a tiny 
expansion of the structure, reducing the Mn-d (e) band 
width, will move the Fermi energy into the minority spin 
gap. In fact we have predicted |Q that for a lattice 
constant larger than ao=5.8A zincblende MnAs will be a 
half- metal. This result, confirmed by other calculations 
p8| , p9| , |6l| , is very attractive since half-metallic systems 
are sources of fully spin polarized currents. We therefore 
investigate if there are conditions under which such a 
zincblende phase of MnAs can be grown. 



A. MnAs: Zincblende vs NiAs-type structure 

The main obstacle to the growth of zincblende MnAs 
is its instability with respect to the NiAs structure. 
The NiAs-type structure is a hexagonal structure (space 
group PGa/mmc) with [6]-coordinated Mn. We have 
studied the relative stability of the NiAs-type and the 
zincblende structures by comparing the total energy per 
MnAs pair as function of the unit cell volume |6^. The 
results are presented in figure ^. It is clear that the 
NiAs-type structure has a much lower total energy and 
also a denser lattice. Therefore it is the stable structure 
at all thermodynamically accessible pressures. However 
we note that if the lattice is forced to expand there is a 
crossover between the two structures, with the zincblende 
being favorable for volumes larger than '^47 A'^ per MnAs 
pair. This corresponds to a lattice constant for the 
zincblende structure of ao=5.8 A, the same at which the 
transition to half-metal occurs. Therefore there is the 
hope that zincblende MnAs could be grown if "negative 
pressure" were applied, for example if it were grown on 
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FIG. 2: Total energy as a function of the MnAs pair volume 
for the NiAs-type (squares) and the zincblende (circles) struc- 
ture. Note the large stability of the NiAs-type structure over 
a very broad volume range. In the inset we present the same 
quantity for the NiAs-type structure as a function of c/a. 



a substrate with large lattice mismatch. Unfortunately 
this is also very challenging. 

In the insert of figure ^ we show the total energy as a 
function of the ratio between the two lattice constants of 
the NiAs-type structure (c/a ratio) at constant volume 
(the experimental volume). The figure shows that the 
NiAs-type structure can accommodate large structural 
distortions without corresponding large energy costs (the 
total energy changes of only about 20 meV when going 
from c/a = 1.7 to c/a — 1.4). This means that at equi- 
librium it is energetically more favorable for the system 
to distort the cell, instead of increasing the volume and 
inducing a NiAs-type to zincblende transition. In conclu- 
sion DPT calculations have shown that zincblende MnAs 
would indeed have desirable properties but that it is un- 
stable and will always tends to form the less attractive 
NiAs phase. Results analogous to those presented here 
are obtained for MnBi MnSb and MnP ||||, suggest- 
ing that the low solubility limit of Mn is a characteristic 
of all the III-V's. 



B. The effect of Spin-orbit coupling 

The commonly used p-d interaction Hamiltonian for di- 
lute magnetic semiconductors is of the Kondo form given 
in equation (|^). The strength of the interaction is gov- 
erned by the exchange constant iVo/3. In this Section we 
will address the question of how A^o/? is affected by the 
spin-orbit coupling. 

Despite the formal similarity of equation (|^) to a 
Heisenberg exchange interaction the "p-d interaction in 
DMS materials does not originate from the Coulomb in- 
teraction but arises as a result of the hybridization be- 
tween p and d derived bands in the crystal. The ap- 
propriate Hamiltonian to describe transition metal im- 



6 



purities in a host crystal is the well-known Anderson 
Hamiltonian. However, Schrieffer and Wolff have shown 
that the Anderson Hamiltonian can be transformed into 
a Kondo-like form | |6^ , containing a term similar to that 
of equation ([|) . This transformation relates the effective 
exchange integral which will be negative in general, 
to the matrix elements of the interaction potential be- 
tween the bands of the crystal. We will discuss the results 
of this section using the following model Hamiltonian for 
the valence band states: 



H = Ho + H 



pd 
X 



H,. 



(7) 



Here is the crystal Hamiltonian without p-d and with- 
out spin-orbit interaction. is the p-d interaction 
given by equation (|l|) and Hsoc is the spin-orbit interac- 
tion. 

Naturally, if the potential changes, the hybridization 
between the bands will be affected and thus the effective 
exchange integral will change accordingly. Note that if 
/3 were a "real" Coulomb exchange integral it would not 
be directly affected by a change in the potential. Only a 
very small indirect effect due to the change of the self- 
consistent charge distribution would be expected in that 
case. 

In the (Ga,Mn)As crystal the As anions introduce a 
substantial spin-orbit interaction which is about one or- 
der of magnitude smaller than the exchange constants de- 
termined within scalar-relativistic calculations. In state 
of the art scalar-relativistic density functional calcula- 
tions spin-orbit coupling is taken into account only ef- 
fectively by using averaged potentials, which conserve 
all non-relativistic symmetry properties of the electronic 
states in the crystal. However in practice the symmetry 
relations change drastically when spin-orbit coupling is 
considered for a spin-polarized system. 

There are two questions to be answered. First, is 
it still possible to transform a hybridization interaction 
into the Kondo-form when the potential becomes spin- 
dependent, as is the case when spin-orbit coupling is 
included explicitly? Secondly, if the Kondo-form still 
holds, how much will the exchange constant change due 
to the change in the potential? In order to answer these 
questions we performed density functional calculations 
based on fuUy-relativistic pseudopotentials |^ for GaAs 
and MnAs, using the code Spinor. Naturally the ob- 
tained density functional results are fully self-consistent 
and do not lend themselves easily to discuss individual 
terms of the model Hamiltonian introduced in equation 
(0) separately. We therefore proceed in the opposite di- 
rection and check whether the density functional results 
can be interpreted by the model Hamiltonian. 

We assume that of equation has been solved 
for the valence band states at the Brillouin zone center. 
For scalar-relativistic band structures of zincblende semi- 
conductors like GaAs the valence band top is six times 
degenerate. This is also the case for the magnetic system 
MnAs as long as is turned off. However, when 



is turned on the p-like states at the top of the valence 
band will split according to their spin orientation into 
two groups of triply degenerate states, separated by the 
energy Nq(3{Sz), following equation (1) (see also the dis- 
cussion of Section IV B). Here {Sz) is the z-component 
of the spin-polarization per unit-cell. This is the result 
found by the scalar-relativistic approach. Finally, treat- 
ing Hsoc as a perturbation, all remaining degeneracies 
are lifted by first order energy corrections as is shown in 
figure H Note that the second order corrections are only 
of order 1 meV for typical values of A, the spin-orbit 
splitting of the valence band at F, and B = ^Nol3{Sz)- 
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FIG. 3: Exchange and spin-orbit splitting of the valence band 
top in MnAs. The split of the center of the band is 6B where 

Now let's turn to the density functional results. Table 
I lists the the scalar- and fully relativistic valence band 
edge energies for MnAs. The DFT calculations confirm 
the model discussed above in that all scalar-relativistic 
degeneracies are lifted when spin-orbit coupling is in- 
cluded explicitly. We also find that, in both groups, 
one state lies approximately at the same energy as the 
scalar-relativistic degeneracy, as predicted by the model. 
However, the energy difference between the three minor- 
ity states is approximately 50 meV, whereas the same 
number for the majority states is found to be approxi- 
mately 100 meV. According to figure ^ the model pre- 
dicts a splitting of |A for minority and majority states 
alike, which is approximately 80 meV, using the GaAs 
spin-orbit splitting for A. Hence the model and density 
functional results show significant discrepancies of about 
30% which would indicate that the model Hamiltonian of 
equation (0), and therefore the Kondo-form of the p-d in- 
teraction, is questionable for MnAs when spin-orbit cou- 
pling is taken into account explicitly. However, the de- 
viations are about two orders of magnitude smaller than 
the splitting between the majority and minority states. 
Hence the corrections due to spin-orbit coupling affect 
the exchange constant Nq/J for MnAs only in the order 
of 1%. 

The model Hamiltonian given in equation (|^) thus re- 
mains a good approximation for the high concentration 
limit as was shown for MnAs. However, there are two 
important issues when the Mn concentration is lowered. 
First, as the splitting between the minority and major- 
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Minority ( 


eV) Majority (eV) 


without Hsoc 


-2.405 


2.405 


with Hsoc 


-2.355 


2.499 


with Haoc 


-2.408 


2.399 


with Hsoc 


-2.458 


2.308 



TABLE I: Top of the valence band levels for MnAs at the Bril- 
louin zone center determined by density functional method 
with and without spin-orbit coupling. 

ity valence band top becomes smaller the observed de- 
viations between the model and the density functional 
results will become increasingly important. There might 
be a point reached where the Kondo form of the p-d in- 
teraction given in equation (|l|) breaks down. This argu- 
ment, however, is based on a mean field extrapolation 
into the low density regime. Results presented in section 
IV of this paper indicate that a mean field approximation 
might be invalid for the p-d interaction in dilute magnetic 
semiconductors altogether. 

Second, as the Mn concentration decreases the band 
structure, especially around the top of the valence 
band, becomes strongly dominated by spin-orbit effects. 
Hence, scalar-relativistic results in the low density limit 
need to be interpreted very carefully, especially when 
spin-dynamics is considered. 

IV. Gai^Mn^As IN THE DILUTED LIMIT 

The first principle study of Gai-ajMn^; in the low di- 
lution limit is a formidable theoretical challenge. This 
is due to the large number of atoms that one should in- 
clude in the unit cell in order to reproduce the exper- 
imental Mn concentrations. However several numerical 
implementations of DFT capable of dealing with a large 
number of atoms are now available. These are generally 
based on the use of pseudopotentials and on local- 
ized basis sets |Q . Most of the results we will present in 
the following sections are obtained with the code siesta 
|52| , which combines both of these features. 

A. Electronic and magnetic properties 

We start by analyzing the calculated density of states 
(DOS) of a 64-atom GaAs unit cell containing one 
MuGa substitution (figure |6^. This corresponds to 
a;=0. 03125, an experimentally accessible concentration. 
From the figure it is clear that Gai_a;Mnj;As has the elec- 
tronic structure of a half-metal. This result is largely con- 
firmed by all density functional calculations to date both 
using the LDA g, |9[ |7|] and the GGA |7l| approxima- 
tion. It is also the same structure found for Ini^x^T^x^^ 
[ |72| . If we project the density of states onto the differ- 
ent orbital components (partial density of states, PDOS) 
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FIG. 4: Partial density of states for Gai-ajMn^^As for x = 
0.03125 (one Mnca in a 64 atom GaAs cell): (a) majority 
spin, (b) minority spin. The vertical line denotes the position 
of the Fermi energy. 



some additional features are revealed. The majority band 
exhibits two broad peaks between -4 eV and -1 eV be- 
low the Fermi energy with strong Mn-d e and t2 compo- 
nent respectively. In addition there is a rather narrow t2 
peak at the Fermi energy. In contrast, the minority band 
has almost no d-character below Ep but instead has two 
sharp e and t2 peaks around 1 eV above Ep . The different 
peak widths reflect the different degrees of hybridization 
of the Mn-d sub-bands with the GaAs bands. 

The magnetic moment of the unit cell is 4 /ib, and 
remains the same up to concentrations of the order of 
x=0.5 | |7l[ . An integer number for the magnetic moment 
is consistent with the half-metallicity seen in the DOS. 
At this point it is therefore very relevant to discuss the 
atomic configuration of Mn in GaAs. The band structure 
around the F point for a;=0. 03125 Gai_j;Mn;:cAs (1 Mn 
ion in a cubic 64 atom GaAs cell) is presented in figure 
^, along with the orbital resolved DOS at the F point. 
We consider the band structure only along the direction 
(iiL,o,0) ^ (0,0,0) ^ (liL, liL, liL) with CO the unit 
vector of the cubic cell. We mdicate these two directions 
respectively as X and M. 

From the picture it is very clear that the Fermi energy 
cuts through the top of the valence band for the majority 
spin. This, in addition to the fact that all the Mn-d states 
in the majority spin band are occupied, suggests that 
Mn in GaAs is incorporated as Mn^"*" and that there is a 
polarized hole, which is antiferromagnetically coupled to 
the Mn. The presence of a hole is revealed by an accurate 
analysis of the integrated DOS ||6^ , and a signature of the 
antifcrromagnetic coupling is the fact that the induced 
magnetic moment at the As sites is antiparallel to that 
of the Mn ions ||, 0. 

It is important to stress that the results from DFT are 
not as ideal as the atomic-like picture discussed above. 
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FIG. 5: Band structure and orbital resolved DOS at the T 
point for Gai-iMn^^As with a;=0. 03125 (1 Mn ion in a cubic 
64 atom GaAs cell): (a) majority band, (b) minority. The 
horizontal line denotes the position of the Fermi energy. 



which predicts S = 5/2 for the Mn. MuUiken population 
analysis ^ |7^ shows that the population of the Mn 
d orbital is 4.7 and 0.7 electronic charges respectively 
for the majority and minority states. Considering the 
fact that the overlap population is of the order of 0.7, 
the raw data are compatible with both Mn df^ and d^. 
Moreover we must stress that at the top of the valence 
band there is a quite large hybridization between Mn d 
and As p states. All these factors contribute to reduce the 
magnetic moment from 5 {S = 5/2) expected from the 
atomic-like picture to 4 /^b ■ A magnetic moment of 4 /ie 
appears to be smaller than that found experimentally 
4.4/xb) (T^, although the agreement can be restored 
by considering partial hole compensation as we will show 
in the following sessions. 

The behavior of Mn in GaAs seems is similar to that 
of Mn substituting the cation sites in other III-V, and in 
IT VI and group IV semiconductors. Schulthess and But- 
ler have calculated the electronic structure of Mn in 
GaAs, Ge, ZnSe and ZnO. The main results are that i) 
the number of minority electrons is not changed by the 
Mn impurity and ii) each Mn impurity adds five addi- 
tional majority states to the valence band. This leads to 
a magnetic moment of 3 /^b, 4 /iB and 5 /ib, respectively 
for Mn in Ge, GaAs, and both ZnSe and ZnO. 

It is difficult to extract the localization properties of 
the holes introduced by the Mn ions from DFT calcula- 
tions. In particular we are not able to conclusively estab- 
lish whether or not the holes are bound to Mn^+ forming 
a neutral {3d^ + h) complex |19 . On the one hand plots 



of the charge density obtained from states within 0.5 eV 
around the Fermi energy reveal that most of the charge 
is concentrated around the Mn sites Q. This seems to 
suggest localization of the hole around the Mn ion. On 
the other hand accurate valence band fitting |Q reveals 
an effective mass quite similar to that of GaAs. This is 



of course an indication of delocalization. We believe that 
this point needs further investigation. A key element to 
determining the localization properties of the holes is a 
knowledge of the exchange constant iVo/3, which we dis- 
cuss in the next Section. 



B. The exchange coupling 

We calculate the exchange constant by evaluating the 
spin-splitting of the conduction and valence bands. This 
mimics a typical magneto-optical experiment The 
main idea is that in the mean field theory based on the 
Hamiltonian of equation (|^) the spin-splitting of the va- 
lence band depends linearly on both the exchange con- 
stant NqP and the Mn concentration x |6^. The same 
argument holds for the spin-splitting of the conduction 
band, which is regulated by a similar Hamiltonian with 
exchange constant Naa. Therefore the exchange con- 
stants can be directly computed from the conduction 
band-edge (valence band-edge) spin-splittings Ai?'^ = 



El 



El ~ {AE" 

Noa = AEyx{S) 



i?P as follows 



NoP^AE^xiS), (8) 



where (S) is half of the computed magnetization per Mn 
ion. 

Recalling the fact that (Ga,Mn)As has a direct gap 
at the r point, we calculate the band structure of 
Gai_a,Mn2.As supercells (see figure ^ around the F point 
for different Mn concentrations, and extract the exchange 
constants by using the equations (||). Our results are in 
table ||. 

X AE" (eV) Ag" (eV) Npa (eV) NqP (eV) 



0.06250 


0.0339 


-0.6839 


0.272 


-5.48 


0.04166 


0.0248 


-0.5458 


0.298 


-6.54 


0.03125 


0.0105 


-0.4472 


0.168 


-7.34 


0.02084 


0.0099 


-0.3442 


0.234 


-8.16 



TABLE II: Conduction AE'' and valence AE" band-edge 
spin-splitting, and exchange constants as a function of the 
Mn concentration x for Gai-xMna;As. 

The behaviors of the valence and conduction bands are 
remarkably different. For the conduction band, although 
the spin splitting shows large fluctuations with x, there is 
no systematic variation with the Mn concentration. This 
indicates that the mean field approximation that led to 
equation ^ is appropriate and one can conclude that the 
exchange coupling between electrons in the conduction 
band and the Mn is ferromagnetic with an exchange con- 
stant Noa ~0.2 eV. This is expected since the coupling 
in this case is direct (Coulombic s-d coupling). Note also 
that the value of the exchange constant Na is very close 
to that usually found in II- VI semiconductors ||7^ . 

In contrast the valence band shows strong deviation 
from the mean field expression (ra) , since the valence band 
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spin-splitting does not vary linearly with x. Turning the 
argument around, -/Vo/3 increases with decreasing Mn con- 
centration, a behavior already well known to occur in 
Cdi_2;Mn^S ||7|, m This suggests that the mean 
field approximation leading to the equations (H) is not 
appropriate for the valence band of (Ga,Mn)As. 

A breakdown of the mean field model occurs when the 
potential introduced by the Mn ions is comparable with 
the relevant band-width. We have calculated Q the cor- 
rections to the mean field model by using a free electron 
model with magnetic impurities described by square po- 
tentials. The calculation is based on the theory of Bcnoit 
a la Guillaume, Scalbert and Dietl, who computed the 
energy within the Wigner-Seitz approach ||80| . 

The main result of the corrected theory is that the 
mean field approximation tends to underestimate the ex- 
change coupling for low dilutions, as observed in our 
LSDA calculations. Our estimation of the exchange con- 
stant gives a value in the range of —4.9 eV < NqP < 
—4.4 eV, which is very large if compared with the values 
quoted by experiments. Moreover it is very important 
to note that within this model the valence holes appear 
to be nearly bound to the Mn ions. This adds a further 
indication of the existence of the {Scf' + h) complex, at 
least in the low dilution limit. 

Finally we want to point out the fact that both the 
exchange constant and the valence band spin-splitting 
are much larger than that found in typical experiments. 
There may be several reasons for this disagreement. 
First, there is strong experimental evidence in the ab- 
sence of saturation in the M — H curves at large mag- 
netic fields and in recent x-ray magnetic dichroism 
measurements , that not all Mn ions contribute to the 
ferromagnetism. Of course an overestimation of x leads 
to an underestimation of Nq/]. Secondly the well known 
lack of accuracy of the LDA to describe strongly localized 
charges may result in an over-estimation of the p-d cou- 
pling |8|]. This of course leads to a larger NqP constant. 
However, we have shown that LDA does not strongly 
overbind the NiAs-type MnAs, and that the structural 
properties of zincblende MnAs are very similar if calcu- 
lated with LDA or LDA+U. This suggests that the error 
in the determination of the p-d coupling is not dramatic 
within LDA. Therefore we do believe that the two main 
conclusions from our LDA calculations, namely that the 
exchange constant is large enough for the mean field ap- 
proximation to breakdown, and that the holes are nearly 
bound, are indeed reliable. It is interesting to remark 
that in a recent paper Chattopadhyay, Das Sarma 
and Millis found a Tc for (Ga,Mn)As in very good agree- 
ment with the experiments, by using dynamical mean 
field theory and our value for the exchange constant. 
Moreover the dynamical mean field theory estimation of 
the critical exchange constant needed for the formation 
of an impurity band, agrees very well with our value [ |68[ . 
This confirms the reliability of our analysis in the exper- 
imentally relevant regime. 



V. THE IMPORTANCE OF INTRINSIC 
DEFECTS 

So far we have always considered the ideal case in 
which only Mn ions are introduced in GaAs. If we assume 
the nominal valences for Mn (Mn^+) and Ga (Ga^+), 
we conclude that Mn acts as single acceptor in GaAs. 
Therefore an equal concentration of Mn ions and holes 
is expected. In contrast in the actual samples the hole 
concentration is much smaller than that of Mn Q and 
some compensation mechanism occurs. As suggested in 
the first experimental works, the presence of As antisites, 
AsGa, usually quite abundant in low-temperature GaAs 
Ql, is likely responsible for the compensation (Asca in 
GaAs is a double donor). Here we summarize the effect 
of AsGa on the magnetic properties of (Ga,Mn)As. The 
detailed results, obtained with SIESTA, can be found in 
references [|| |6[ l^. 



A. Compensation due to Asca: local effects 

We consider explicitly the effects of the inclusion of 
AsGa in (Ga,Mn)As at different dilutions, and study how 
the chemical environment modifies the magnetic interac- 
tion between the Mn ions. We construct 64 (cubic) and 
32 (rectangular) atom GaAs cells in which we include two 
Mn ions (leading to Mn concentrations of respectively 
x=0.0625 and a;=0.125) and a variable number of Asoa 
antisites. The parameter which quantifies the strength 
of the exchange coupling is the energy difference ApA 
between the total energies of the antiferromagnetic Eaf 
and ferromagnetic EpM configurations of the supercell. 
These are obtained by fixing the spin direction at the 
beginning of the self-consistent calculation. 

We further look at two possible spatial configurations 
of the Mn ions in the unit cell: 1) separated, when the 
Mn ions occupy positions as far apart as possible (i.e. the 
corner and the middle of the cubic cell), 2) close, when 
the Mn atoms occupy two corners of a tetrahedron and 
are coordinated through a single As ion (see figure 

ApA, and the magnetization per Mn ion, Mmh, for 
the separated arrangement are presented in figure || as a 
function of the number of Asca antisites. The magne- 
tization per Mn ion is defined to be half of the magne- 
tization of the cell, calculated in the FM aligned phase. 
First we note that the ferromagnetic coupling is strongly 
weakened by Asoa antisite doping. This is expected ac- 
cording to the picture of hole-mediated ferromagnetism: 
AsGa antisites contribute electrons into the system and 
therefore compensate the holes. We also note that in 
the case of no antisites, ApA is much larger in the case 
of large Mn concentrations (we recall that according to 
equation (|^) Tc scales linearly with x). The figure also 
suggests that the compensation mechanism does not fol- 
low the nominal atomic valence, since a single Asca anti- 
site per cell is not sufficient to destroy the ferromagnetic 
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FIG. 6: (a) Energy difference between AF and FM align- 
ments, ApA, and (b) magnetization per Mn ion, Mmu, as a 
function of the number of Asca antisites in the cell: sepa- 
rated configuration. The horizontal line denotes the division 
between FM and AF alignment. 



coupling. Above compensation (more that one Asca for 
two Mn ions) antiferromagnetic coupling is obtained for 
large Mn concentration, while the system stays ferromag- 
netic at low concentration, although in both cases Afa is 
rather small (|Afa| < 20 meV). This is consistent with 
the onset of antiferromagnetic super-exchange coupling 
p8[ , the mechanism which is believed to be responsible 
for the magnetic order in the II- VI DMS |76) , at compen- 
sation. Super- exchange is a short range interaction and 
therefore is less important in the low concentration limit 
where the Mn ions are well separated. 

These results agree qualitatively with those obtained 
by Akai ^ for (In,Mn)As, using a KKR-CPA-LDA 
(Korringa-Kohn-Rostoker Coherent Potential Approxi- 
mation and Local Density Approximation) method 
Akai interpreted his data as a competition between ferro- 
magnetic double-exchange and antiferromagnetic super- 
exchange. 

In figure ^ we see that Mmh increases with the 
AsGa concentration and than saturates to a value around 
4.5 fiB- This can be easily explained by remember- 
ing that the top of the majority spin valence band of 
(Ga,Mn)As has some Mn d component due to hybridiza- 
tion (see figure H). Asca-doping moves the Fermi energy 
towards the conduction band, filling the valence band 
completely. This enhances the magnetic moment of the 
unit cell. The saturation is due to the fact that the next 
Mn d states available above the valence band are at the 
edge of the conduction band in the minority spin band. 
However Asca doping pins the Fermi energy at mid gap 
and these states can not be filled. Note that a magnetic 
moment per Mn of 4.5 /is is in good agreement with 
x-ray circular magnetic dichroism measurements | |75| . 

A better insight into the mechanism giving rise to the 
ferromagnetic order is given by the results for the close 
configuration. In this case the local microscopic configu- 
ration is crucial for the magnetic properties, therefore we 
consider three different situations (see figure ^): a) the 
antisites are far from the Ga2Mn2Asi complex, b) one 
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FIG. 7: Energy difference between the AF and FM alignments 
ApA (left), and magnetization per Mn ion Mmii (right) as a 
function of the number of Asca antisites in the cell: Close Mn 
arrangement. The symbols •, * and o represent arrangements 
(a), (b) and (c) respectively. 



antisite occupies a tetrahedral site (GaiMn2As2), and c) 
two antisites occupy the tetrahedral sites (Mn2As3). The 
results for a;=0.0625 are presented in figure |^. 

In case (a) both the ApA and Mun curves look very 
similar to those found for the separated arrangement (fig- 
ure ^ . The main difference is a strong antiferromagnetic 
couphng above compensation (n Asca >!)■ RecalHng 
that Mn in GaAs assumes the d^ configuration (rf band 
half- filled), and that the Mn d shell is antiferromagneti- 
cally coupled with the intermediate As atom (through the 
p-d interaction), we conclude that super-exchange cou- 
pling stabilizes the AF phase at and above compensation 
in the Ga2Mn2Asi complex. Note that the same con- 
clusions have been drawn by Park et al. for Gei-xMji^ 
pl| , although in that case the antiferromagnetic coupling 
seems to occurs also if holes are present. 

Cases (b) and (c) present several interesting features. 
The most remarkable is that the ferromagnetic alignment 
is stable and almost insensitive to the total Asca con- 
centration. This suggests that the dominant interaction 
in the GaiMn2As2 and Mn2As3 complexes is completely 
local. Once again MuUiken population analyses help in 
understanding the electronic configuration of these com- 
plexes. In table III we present the MuUiken orbital popu- 



n AsGa Type Mn-d| 


Mn-d| 


As-pi 


As-pi 


As-p 


As 


2 a 


4.74 


0.70 


1.55 


1.63 


3.18 


4.92 


2 b 


4.74 


0.71 


1.50 


1.67 


3.17 


4.95 


2 c 


4.72 


0.75 


1.44 


1.71 


3.15 


4.97 


3 a 


4.76 


0.68 


1.57 


1.63 


3.20 


4.93 


3 b 


4.74 


0.74 


1.55 


1.65 


3.20 


4.97 


3 c 


4.73 


0.75 


1.49 


1.69 


3.17 


4.99 



TABLE III: MuUiken atomic and orbital populations for the 
Mn ions and the intermediate As atom of the complexes of 
figure ^. The Mn concentration is a;=0.0625. The symbols t 
and J, correspond to majority and minority spin respectively. 
The populations are in units of the electronic charge |e|. 
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lation for the two Mn ions and the intermediate As ion of 
the complexes of Fig.^. To understand the table, imagine 
the situation in which two Asca antisites are first located 
far from the Mn-As-Mn complex, and then each antisite 
is moved in turn to one of the two other corners of the 
tetrahedron. This corresponds to moving vertically (top 
to bottom) down the first half of the table. We notice 
that: i) the charge on the middle As atom increases, ii) 
the spin-polarization of the p shell of the middle As atom 
increases, and iii) the total population of the p-shell of 
the middle As atom decreases. And most importantly 
the magnetic coupling changes from antiferromagnetic to 
ferromagnetic. 

It is well-known that in an antiferromagnetic crystal 
the presence of a bound carrier (electron or hole) which 
is Zener coupled to the local spins always induces a dis- 
tortion in the antiferromagnetic lattice |Q. We there- 
fore propose that the observed transition from antifer- 
romagnetic in Ga2Mn2Asi to ferromagnetic coupling in 
GaiMn2As2 and Mn2As3 results from the onset of ferro- 
magnetic double-exchange coupling mediated by a bound 
Zener carrier. It is important to note that the potential 
necessary to bound an extra charge in the vicinity of the 
Mn ions is provided by the As antisites. We also empha- 
size that if Ga2Mn2Asi or Mn2As3 complexes are present 
in actual samples, then the Mn ions in the complexes will 
not contribute to the overall ferromagnetic alignment, 
being magnetically "locked" by the local environment. 
This is a possible explanation of the fact that in actual 
samples only a fraction of the Mn ions contribute to the 
ferromagnetism 0. Similar conclusions based on 
DFT calculations have been drawn for Gei-xMux |ll|. 





FIG. 8: (a) Asca and (b) Asi-Voa pair obtained by moving 
AsGa along the (111) direction. 



as a function of the displacement ^(iii) of Asca along 
(111). For comparison we also present the same curve for 
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FIG. 9: Total energy (left-hand side scale) and magnetization 
(right-hand side scale) for (a) (Ga,Mn)As and (b) GaAs as a 
function of the displacement of an Asca antisite along (111). 
The energy of Z(iii)=0 is set to eV. 



B. Defect manipulation 

In a recent paper [ p6[ we have shown that intrinsic de- 
fect manipulation can provide a valuable way of tuning 
the carrier concentration and hence enhancing the Tc in 
(Ga,Mn)As, without changing the Mn concentration nor 
the microscopic configuration of the Mn ions. Here we 
summarize the basic ideas. We first recall that isolated 
AsGa in GaAs are responsible for the photoquenchable 
EL2 defect In fact upon illumination Asca under- 

goes a structural transition to an As interstitial-Ga va- 
cancy (Asi-Vca) pair, which is obtained by moving Asca 
along the (111) direction (see figure This complex is 
metastable since Asoa can be regenerated by heating. It 
is crucial to observe that the Asj-Vca pair is not electron- 
ically active in GaAs, since its only state in the bandgap 
is completely filled. In what follows we show that this 
metastable complex is present, can be obtained by illumi- 
nation and is also electronically inactive in (Ga,Mn)As. 

In figure |^ we present the total energy and the mag- 
netization of a 64 atom unit cell containing one Mn ion 
(x=0. 03125) and one Asca (note that in this case the 
system is n-doped, with one excess electron per unit cell) 



GaAs, which is in very good agreement with previously 
published results It is clear that also in (Ga,Mn)As 
the Asi-Vca defect is metastable and that the energy 
barrier for the thermal regeneration is 0.45 eV. This, ac- 
cording to kinetic calculations |9l| gives a regeneration 
temperature of about 100 K. 

Having established that the AsrVca defect is 
metastable let us now prove that it can be generated also 
in (Ga,Mn)As. The mechanism has been explained by 
Scheffier et al. ^ , and can be summarized as follows. 
It is known that the excited a\t\ electronic configuration 
of a tetrahedral substitutional double donor induces lat- 
tice distortion. Such distortion occurs is Asca, which 
possesses a doubly occupied donor level a at midgap and 
an empty resonant state with t2 symmetry close to the 
conduction band edge. The distortion is initiated be- 
cause the many-electron wave function of the a\t\ con- 
figuration is orbitally degenerate therefore the system is 
Jahn- Teller unstable. Jahn- Teller distortion splits the 
t2 state into a lower a state (half-filled) and a higher e 
state (empty). Therefore an optical excitation of the a\t^ 
ground-state to the a\t\ will initiate a distortion. It has 
been demonstrated that in GaAs the total-energy curve 
as a function of the displacement of Asca along (111) for 
the a\t\ configuration has a minimum for ^(iii) ~0.3A 
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pl| . Then the system has some probabihty of relaxing 
onto the Asi-Vca side of the total-energy curve of figure 
^, ending up considerably far from the antisite position, 
creating the Asi-Vca pair. 

The same mechanism holds for (Ga,Mn)As since the 
only difference with respect to GaAs is the spin-splitting 
of the bands. This splitting however does not move the 
a state of Asoa into the valence band nor the ^2 state 
into the conduction band. With these considerations it is 
clear that the As antisite in (Ga,Mn)As presents the same 
features as that in GaAs and, therefore, the mechanism 
described above is still applicable. 

The effects of this transition on the ferromagnetic 
properties are quite dramatic. We calculate the energy 
difference between the ferromagnetic and antiferromag- 
netic alignment, ApA, and we observe an increase from 
53 meV to 124 meV when a single AsQa is transformed 
into a Asi-Vca pair. Hence the ferromagnetic order is 
strengthened. The reason for this enhancement of the 
ferromagnetic coupling is that most of the hole compen- 
sation introduced by Asca is lifted by the creation of the 
Asi-Vca pair- The band structures of (Ga,Mn)As con- 
taining AsGa and Asj-Vca pair shown in figures |l^and|ll| 
explain this point. It is clear that in the case of Asca the 
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FIG. 10: Band structure for Gai-ajMn^^As with intrinsic de- 
fects: AsGa- These are the bands for a cubic 64 atom supercell 
(2x2x2 zincblende cubic cells) containing one Mn ion and one 
AsGa- On the right-hand side the majority spin and on the 
left the minority. The horizontal line denotes the position of 
the Fermi energy. 

Fermi energy cuts through the Asca impurity band, leav- 
ing the GaAs valence band completely filled (note that 
the donor level is only marginally spin-split, although it 
is hybridized with the top of the valence band). This 
means that the compensation is complete (note that in 
the figure we have one Mn ion and one Asca in the cell, 
therefore we correctly expect an excess of one electron 
in the impurity band). In contrast in the case of the 
Asi-Vca pair there are still holes in the system, since the 
Fermi energy cuts through the top of the valence band 
and the impurity level. It is also important to note that 
the impurity level has small hybridization with the ma- 




FIG. 11: Band structure for Gai-ajMn^^As with intrinsic de- 
fects: Asi-Voa pair. These are the bands for a cubic 64 atom 
supercell (2x2x2 zincblende cubic cells) containing one Mn 
ion and one Asi-Vca pair. On the right-hand side the ma- 
jority spin and on the left the minority. The horizontal line 
denotes the position of the Fermi energy. 



jority valence band of (Ga,Mn)As. This means that in 
this case there are holes in the system with a dispersion 
very close to that of the defect free case. For this reason 
the presence of As^-Vca pairs does not alter the magnetic 
properties compared with the defect-free alloy (note that 
ApA for this situation is 124 meV, very similar to the 
defect free value of 159 meV). 

In conclusion, the optically induced transition of 
AsGa to a Asi-Vca pair reduces the hole compensation, 
strengthening the ferromagnetic coupling. This, in prin- 
ciple, allows tuning of the hole concentration without act- 
ing on the chemical composition of a sample. However, 
since the antisites regenerate for temperatures of the or- 
der of lOOK, the mechanism cannot be used to obtain 
very high TcS. 



VI. DILUTED FERROMAGNETIC 
HETEROSTRUCTURES: MAGNETIC AND 
TRANSPORT PROPERTIES 

We have seen in the previous Sections that formation 
of the NiAs-type phase of MnAs is easier than forma- 
tion of the zincblende phase, which in turn is very de- 
sirable due to its half-metallic band structure. It has 
been demonstrated |Q, however, that zincblende MnAs 
can be grown up to half of a monolayer in GaAs/MnAs 
superlattices. These structures are called digital ferro- 
magnetic heterostructures (DFH). Since our earlier work 
has shown a strong dependence of the magnetic proper- 
ties on the local arrangement of the Mn ions, we expect 
that the properties of DFHs would be very different from 
their random alloy counterparts. 

These are the main experimental findings p2| : first, Tc 
decays with increasing GaAs layer thickness separating 
the MnAs sub-monolayers, and saturates for thicknesses 
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larger than ~50 GaAs monolayers. The saturation is 
unexpected according to the mean field model for three 
dimensional systems, since the total Mn concentration in 
the sample decreases with the increase of the GaAs thick- 
ness. This separation dependence suggests that DFH's 
behave like planar systems. 

Secondly, Hall measurements in the direction paral- 
lel to the MnAs planes show an anomalous Hall effect 
for undoped samples, which disappears upon Be-doping 
|p3| , 0. Large Shubnikov de Haas oscillations are found 
in doped samples, although surprisingly the charge den- 
sities extracted from the Hall coefficient and from the 
Shubnikov de Haas oscillations are different. This sug- 
gests that two different carrier types could be present in 
the system. 

Density functional theory has answered the following 
questions regarding the physics of DFHs: i) what is the 
real dimensionality of the system? ii) are the carriers 
spin-polarized? iii) what is the carrier distribution in the 
system? 

We have used the code siesta with a DFH su- 
perlattice constructed from N GaAs cubic cells (8 atoms 
in the cell) aligned along the z direction. One Ga plane 
(two atoms) is substituted with Mn and periodic bound- 
ary conditions are applied. This leads to an infinite 
MnAsi/GaAs2Ar-i superlattice, where MnAs zincblende 
monolayers are separated by a 5.65 x A^A thick GaAs 
layer. In figure |l^ we present the calculated band struc- 
ture for the case N = 8 (with a 45. 2 A thick GaAs inter- 
layer), for both the majority and minority spins. It shows 
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FIG. 12: Band structure for a MnAsi/GaAsis superlattice: 
(a) majority and (b) minority spins. Xi and X2 are in the 
MnAs plane and denote directions along respectively the edge 
and the diagonal of the cubic supercell. F — > is the di- 
rection orthogonal to the MnAs plane. The horizontal line 
indicates the Fermi energy. 

clearly a very peculiar half-metallic structure. Although 
the band structure has a gap for the minority spin and 
some band crossing at the Fermi energy for the majority. 



the band dispersion in the MnAs plane {Xi F ^ X2) 
is quite broad while the one perpendicular to the MnAs 
plane {Zi —>■ F) is very narrow (impurity- like band). 
Therefore MnAs / GaAs looks like a two dimensional half 
metal with small hopping between the MnAs planes. 

We also look at the stability of the ferromagnetic phase 
as a function of the separation between Mn planes, by 
calculating ApA- Surprisingly ApA is independent of the 
GaAs thickness for the range of thicknesses investigated 
here (531 meV, 533 meV and 515 meV respectively for 
N—A, N—6 and N=8). This is consistent with the exper- 
imental insensitivity of Tc to the GaAs thickness. This 
first analysis shows that most of the physics of DFH oc- 
curs in the MnAs planes. 

In order to have a better understanding we have also 
performed transport calculations in the ballistic limit. 
One of the advantages of using a localized basis set for the 
DFT calculation is that, at the end of the self-consistent 
procedure a tight-binding Hamiltonian is generated by 
direct numerical integration over a real space grid [ p3[ . 
Then transport properties can be calculated in the bal- 
listic limit by using a well established technique for tight- 
binding Hamiltonians, which is described in reference 
p6| . For this work, we have generalized the technique to 
the case of non-orthogonal tight-binding models with sin- 
gular coupling matrices More details can be found 
in the cited literature. 

In figure |l^ we present our calculated conductance per 
unit area as a function of the position of the Fermi energy 
for a MnAsi/GaAsi5 superlattice for both the current in 
the MnAs plane (GIF) and current perpendicular to the 
Mn plane (GPP) directions, and for both spins. We also 
project the conductance onto the atomic orbital basis set 
in order to determine the orbital character of the elec- 
trons carrying the current fo^. As expected from the 
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FIG. 13: Total and partial conductance per unit area for a 
MnAsi/GaAsi5 superlattice as a function of the position of 
the Fermi energy for majority (|) and minority (J,) spin bands. 
The vertical line denotes the position of the Fermi energy for 
undoped samples. 
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band structure, the ballistic current is 100% spin polar- 
ized, with no current for the minority spin band. In the 
majority spin band the behavior is very different for the 
CIP and CPP alignment. In the CIP case the conduc- 
tance is quite large and independent of the energy, with 
a significant contribution (roughly 20% of the total con- 
ductance) coming from the Mn d t2 orbitals. In contrast, 
the conductance is very small in the CPP direction with 
orbital contribution almost entirely from the As p states. 
Moreover the conductance at comes from just a few k- 
points around the F point (in the direction orthogonal to 
the transport). These correspond to the states with the 
largest kinetic energy in the direction of the transport. 
This situation is similar to that occurring in tunneling 
junctions and so we describe the transport as tunncling- 
like, meaning that the transport in the CPP direction is 
through hopping between the MnAs planes. 

Finally we investigate the spatial distribution of the 
current. This is given by the charge density distribution 
in real space, p(r), calculated only for those states con- 
tributing to the conductance (see figure |lj) |9^. The 




FIG. 14: Charge density distribution in real space for 
MnAsi /GaAsis (a) calculated only for those states contribut- 
ing to the conductance and energy within 0.3 eV from Ep: 
CIP (b) and CPP (c) configurations. 

figure confirms that the current in the CIP case is dis- 
tributed mainly in a narrow region around the MnAs 
planes, with small spillage outside. In contrast, the CPP 
current is mainly located at the Mn plane with small con- 
tributions from the GaAs layers. This means that carri- 
ers are strongly confined in the MnAs plane and the per- 
pendicular transport is via hopping between the planes. 

DFT has therefore answered all the questions stated 
at the beginning of this Section. The DFHs appear to 
be two dimensional half-metals, with a metallic-like bal- 
listic conductance in the MnAs plane and an hopping 
conductance perpendicular to the MnAs planes. The fer- 
romagnetism is therefore insensitive to the GaAs layer 



thickness, showing no dilution effect. 



VII. OTHER DILUTED MAGNETIC 
SEMICONDUCTORS 

So far we have concentrated primarily on the prop- 
erties of (Ga,Mn)As, which is the most well studied of 
the DMSs. However, in the last year several other DMSs 
have been synthesized, some of them showing remarkably 
high Tc- At the same time several DFT calculations have 
been published both explaining the properties of existing 
materials, or making predictions for new ones. The pre- 
diction of material properties in advance of experiments 
is one of the most important aspects of density functional 
theory. Here we review the principal works in this area. 



A. Gai :iM:iAs with M=V, Cr, Fe 

Doping of GaAs with transition metals other than Mn 
has been investigated theoretically using both the full 
potential linearized augmented plane wave (FLAPW) 
method |^ and the atomic sphere approximation Jzo] . 
Since the first method is computationally more demand- 
ing the author concentrates on the electronic structure of 
hypothetical zincblende MAs with M=V, Cr, Mn, Fe, and 
on the high concentration limit of Gai_a:Ma;As (a;=0.25 
and a;=0.125). The main result is that, while the ground 
state of FeAs is antifcrromagnetic, VAs, CrAs, and MnAs 
at the respective equilibrium lattice constants appear to 
be half-metallic ferromagnets with magnetic moments re- 
spectively of 2 /iB, 3 /iB and 4 fj.^ per formula unit. 

Turning attention to Gai_xM^As it is interesting to 
note that the calculated DOS for Gai_j,V2,As in the large 
concentration limit reveals very small hybridization of 
the V 3d orbitals with the GaAs valence band, suggesting 
that in this compound the magnetic coupling may be 
rather small. In contrast, Gai-xCr^jAs shows a DOS 
very similar to that of Gai_a;Mn^As. This suggests that 
Gai_2;Cra;As may be a potential candidate for a high 
Tc DMS. In addition van Schilfgaarde and Mryasov fit 
their DFT results with a pairwise Heisenberg-like energy 
E = — '^ij JijSi ■ Sj, and found that the values of J for 
Gai_2.Cra;As are larger than those of Gai_xMn2:As. This 
suggests stronger ferromagnetism for Gai-_2:Cra;As than 
for Gai_2;Mna;As. However this prediction is strongly 
affected by the fact that the DFT results strongly deviate 
from the Heisenberg-like form. 

Experimentally both Gai-^^Fcj^As and 
Gai_a;Cra;As ||98| have been synthesized and there 
is a little evidence for ferromagnetic order. However very 
recently zincblende CrAs has been successfully grown on 
GaAs at low temperature , showing a Tc above room 
temperature. Moreover preliminary measurements show 
a magnetic moment per formula unit very close to 3 ug , 
in good agreement with the theoretical predictions [pi]. 
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To our knowledge Gai_a;V2;As has never been grown. 



B. Gai_,Mn,N 

There has been a large effort in synthesizing 
Gai_a;Mn2:N in the last two years, since Dietl et al. pre- 
dicted a Curie temperature higher than any other Mn- 
doped semiconductors for this compound JTzt . Very re- 
cently this prediction has been confirmed experimentally 
pof with Gai_j:Mna;N at 10% Mn concentration show- 
ing a remarkably large Tc (940 K). Unfortunately the 
mechanism for the ferromagnetism in Gai_2:Mna;N is not 
clearly understood due to the lack of a experimental data. 
In particular the nature and the role of the free carriers 
(if present) are not yet clear. 

Sato and Katayama-Yoshida performed DFT calcula- 
tions within the Korringa-Kohn-Rostoker method and 
the coherent potential approximation (KKR-CPA) [|lOO[ . 
They looked at transition metal doping of GaN and 
studied the stability of the ferromagnetic configuration 
against the spin-glass state. The ferromagnetic state ap- 
pears to be stable for (Ga,V)N and (Ga,Cr)N, while for 
(Ga,Fe)N, (Ga,Co)N and (Ga,Ni)N the spin glass state 
has lower total energy at every concentrations studied 
(x=0.05, 0.1, 0.15, 0.20, 0.25). The case of (Ga,Mn)N is 
critical since the ferromagnetic state is stable only at low 
Mn concentration {x <0.15). This is explained in terms 
of competition between double- and super-exchange, with 
the latter being dominant for small Mn-Mn separations. 

The DOS of (Ga,Mn)N looks rather different that of 
(Ga,Mn)As, showing a much stronger d contribution at 
the Fermi energy and smaller p-d hybridization. This 
leads to the formation of a Mn d impurity band i n th e 
band gap of GaN, as calculated by Kronik et al. [101|, 
who also calculated that the (Ga,Mn)N valence band is 
not spin-split. The atomic configuration of Mn is cal- 
culated to be Mn^+ with the d orbitals arranging as d'^ 
[ [loot • This contradict s the experimental results for para- 
magnetic (Ga,Mn)As [102 , which convincingly show a d^ 
{S — 5/2) configuration. The reasons of this discrepancy 
are not clear at the moment. 



C. Zni_:,M^O with M=V, Cr, Mn, Fe, Co, Ni 

The mean field model of Dietl et al. also predicts a 
very large Tc for Mn-doped ZnO |l^, provided that the 
sample is p-doped. Sato and Katayama-Yoshida also in- 
vestigated the stability of the ferromagnetic phase with 
respect to the spin-glass phase for Zni_a;Mj,0 (M=V, 
Cr, Mn, Fe, Co, Ni), again using a KKR-CPA method 



[ 103| , |104| , |105[ . They found a ferromagnetic ground state 
for all the materials except Mn, for which the spin-glass 
configuration has a lower energy. However they also 
showed that in the case of Mn the ferromagnetic configu- 
ration can be obtained by large p-doping (they substitute 



N atoms at the O sites in their calculations). 

Since in the II- VI semiconductors only very small hole 
concentrations can be obtained by doping, they also in- 
vestigated the conditions to have strong ferromagnetism 
in transition-metal-doped ZnO with additional n dop- 
ing (Ga at the Zn sites). The result of this calculation 
is that the presence of electrons stabilizes the ferromag- 
netic order in ZnO doped with Fe, Co and Ni [104[. It 
is important to note that in all these cases the Mn d 
shell is more than half-filled and that the Fermi energy 
is located within a Mn d region of the minority spin 
density of states. These are the conditions for strong 
double-exchange coupling between the transition metal 
ions. The fact that n-doped (Zn,Co)0 [0 and (Zn,Ni)0 
[106[ have been grown, showing TcS above room temper- 
ature, is very encouraging. 

Similar calculations have been carried out for 
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transition- mctal-dopcd ZnS, ZnSe and ZnTe [ 107 
without any additional doping. These show that only 
V- and Cr- doped materials have a ferromagnetic ground 
state, which is due to double-exchange coupling. In this 
case the carriers mediating the double-exchange are holes 
at the top of the valence band, which in turn are strongly 
hybridized with the p states of the group VI element. It 
is worth mentioning that to date no ferromagnetic order 
has been found experimentally in these materials. How- 
ever real samples present strong self-compensation, which 
suppresses the hole-mediated double-exchange mecha- 



D. Mi_^Mn^GeIl2 with M=Cd, Zn and II=P, As 

Cdi_a;Mna;GeP2 was the first room temperature di- 
luted ferromagnetic semiconductors to be grown [|l^. It 
has a body-centered tetragonal cell and the main ad- 
vantage of this structure is that Mn can substitute for 
the II cations, adopting the Mn^+ state ("natural" for 
Mn). Accurate calculations both with GGA and LDA 



have been performed [109, llOj, showing that indeed Mn 
assumes a Mn^+ state. The calculations demonstrate 
a very weak sensitivity of the magnetic and electronic 
structure on both the anion (As or P), and cation (Cd or 
Zn) elements. 

Interestingly both the LDA and GGA results give an 
antiferromagnetic ground state, contradicting the exper- 
imental results. However two aspects must be pointed 
out. First, the samples in the experiments [T^ are pre- 
pared with vacuum deposition on a single crystal surface, 
followed by solid state reaction at high temperature. This 
means that most of the Mn ions are located close to the 
surface. The DFT calculations are performed for perfect 
crystalline bulk systems, whose properties will certainly 
differ from those of the actual samples. Secondly it has 
been demonstrated [110[ that the ferromagnetic phase 
can become stable against the antiferromagnetic one, if 
electrons are introduced in the system (S substituting for 
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P). This seems to suggest that, the presence of free carri- 
ers is also essential for ferromagnetic order in chalcopyrite 
semiconductors. 

Finally, very recently Mahadevan and Zunger calcu- 
lated the formation energies of several intrinsic defects 
of CdGeP2, of Mn impurities in CdGeP2 either a t th e 
Cd and the Ge site, and of complexes of these [111|. 
They found that under Cd, Mn and P rich growth con- 
ditions Mnce impurities naturally form in the crystal. 
These impurities, that in contrast to Mn in GaAs do not 
form clusters, are acceptors in CdGeP2 and align fer- 
romagnetically. In addition it is important to observe 
that the ferromagnetic coupling is also found between 
Mnce and Mncd impurities. These predictions offer very 
good guidelines for growing room-temperature Mn-doped 
CdGeP2. 



VIII. CONCLUSIONS AND FUTURE 
DIRECTIONS 

In this paper we have presented the recent DPT con- 
tributions to the physics of diluted magnetic semicon- 
ductors. Most of the results are for (Ga,Mn)As, have 
been obtained within either the LSDA or the GGA ap- 
proximation and capture most of the relevant physics of 
(Ga,Mn)As. We have also summarize the latest results 
for other DMSs, whose properties are less well studied 
and understood. 

During these studies several disagreements with exper- 
iments have appeared, suggesting that in some cases the 
LSDA/GGA approximations may be not completely ap- 
propriate. It is therefore indispensable to understand the 
limit of the reliability of the LSDA/GGA approximations 
and eventually to correct them. It is well known that the 
local density approximation fails in describing systems 
where electron-electron interaction is strong. This gen- 
erally happens when we try to describe orbitals tightly 
bound to their nuclei. The d orbitals of transition metals 
in transition- metal oxides are an example of this. 

A typical signature of the inappropriateness of the 
LSDA/GGA is that the bands from tightly bound or- 
bitals are pushed towards higher energies. This can cre- 
ate the following problems in DMSs: 1) underestimation 
of the true band occupation, leading to an erroneous mag- 
netic moment, 2) overestimation of p-d interaction, due 
to artificially strong hybridization, 3) overbinding of the 
structure, 4) tendency to metallicity in otherwise semi- 
conducting systems, 5) prediction of the wrong magnetic 
state (ie ferromagnetic instead of antiferromagnetic). All 
these problems can of course strongly affect the predic- 
tions of DPT calculations, and they must be corrected. 

Several schemes have been adopted to go beyond 



the LSDA/GGA approximation, and here we briefly 
list them. In 1981 Perdew and Zunger showed 
that most of the problems of the LDA approximation 
come from its incorrect treatment of coulombic self- 
interaction. They demonstrated that self-interaction cor- 
rections (SIC), namely the subtraction of spurious self- 
interaction from the density functional, greatly improve 
the description of localized stat es. S ince then several im- 
plementations have been made |ll2[ including algorithms 
using pseudopotentials |113|. 



A quite different approach to this problem is repre- 
sented by the LDA+U scheme [|| If. The main idea 
is to combine density functional theory with an Hubbard 
description of those orbitals which suffer strong corre- 
lations. In this case an additional Hubbard-like term 
is added to the energy functional and the consequent 
single particle equations are deduced. The theory gain 
two extra parameters (the Hubbard parameter U and 
the Hund's rule exchange J), which can be fixed from 
experimental values or calculated self-consistently. The 
advantages of the LDA-I-U scheme are the relatively easy 
numerical implementation and the fact that the param- 
eters can be directly compared with those coming from 
simpler models. 

In conclusion, we have shown that DFT can be a very 
powerful tool for describing structural, electronic, mag- 
netic and transport properties of diluted magnetic semi- 
conductors. In particular DFT is very useful for obtain- 
ing parameters not accessible from experiments, and for 
providing accurate descriptions of properties occurring at 
the atomic scale. The prediction of new materials with 
novel properties is possible, although special care should 
be taken in order to keep errors coming from the local 
density approximation under control. 
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